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We investigate in detail the Casimir torque induced by quantum vacuum fluctuations between two 
nanostructured plates. Our calculations are based on the scattering approach and take into account 
the coupling between different modes induced by the shape of the surface which are neglected In any 
sort of proximity approximation or effective medium approach. We then present an experimental 
setup aiming at measuring this torque. 
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Quantum vacuum fluctuations produce Gasimir forces 
when scattered by external boundaries. These bound¬ 
aries may be macroscopic objects such as plates or micro¬ 
spheres, leading to the Gasimir effect [1] or they may be of 
microscopic origin such as atoms or molecules experienc¬ 
ing Van der Waals forces. The dynamics of the Gasimir 
effect between macroscopic boundaries is directly related 
to their geometries, opening a rich physics with a variety 
of extremely interesting theoretical predictions |2]. The 
interplay between geometry and quantum vacuum fluc¬ 
tuations is a topic of high current interest that could lead 
to new tests of Quantum Electro-Dynamics (QED). 

Besides fundamental issues, the interplay has also 
relevant implications in the applied domain. Impor¬ 
tant technological progress has lately been achieved at 
sub-micrometer scales. Nanofabrication techniques of¬ 
fer today the possiblity for designing operational micro- 
and nano-electromechanical systems (MEMS and NEMS) 
[3]. As the dominant interaction between two neutral 
non-magnetic objects in the sub-mirometer range, the 
Gasimir energy is now recognized as a new mean for ac¬ 
tuation, control or functionalization in the actual design 
of MEMS and NEMS |3]. New possibilities of nano¬ 
mechanical control are opened when shaping the inter¬ 
acting objects 1^. 

Due to the dependence of the Gasimir energy to bound¬ 
aries, the source of control is not limited to the normal 
Gasimir effect between planar dielectric or metallic plates 
but has been extended recently by the observation of the 
lateral force between corrugated plates HIT]. This lateral 
force stems from the translational symmetry breaking 
along the plates by the imprint of parallel and periodic 
corrugations on one or both plates. It has been evaluated 
in the past for perfect boundaries using the path integral 
approach ElEl as well as for more realistic dielectric ma¬ 
terials |10l [TT] within the scattering approach |12l Us]. 
The latter allows to compute Gasimir forces between ar¬ 
bitrarily shaped bodies taking into account the material 
properties. The lateral force has been shown to be max¬ 
imized if rectangular structures are imprinted |14| . The 
lateral force has been proposed to be employed in non¬ 


contact rack-and-pinion devices with high potential for 
applications [T3HT7] . 

Recently, exact expressions of the Gasimir force be¬ 
tween dielectric nanostructures have been developed for 
the normal and the lateral regime mnBHii] that com¬ 
pare very well with recent experiments performed with 
deep |22l 133] and shallow corrugations [ 23 ] . 

A different situation of great interest is encountered 
when the parallel and periodic imprints between two 
plates are tilted by a small angle, such as shown in Fig. 
^ In this case a vacuum induced Gasimir torque arises, 
which is the main issue of the present paper. Interest¬ 
ing perspectives for Gasimir torques lie in the design of 
noncontact gears |251 123] which would allow for torque 
transmission between two nanostructured plates avoid¬ 
ing any direct contact between them. 



FIG. 1. Sketch of the configuration giving rise to a Casimir 
torque. 

As the rotational symmetry is broken in this case, the 
Gasimir torque arises because the Gasimir energy now 
depends on the angle 0 between the two directions of the 
corrugations. A preliminary study of this situation has 
already been performed for corrugations small compared 
to the other length-scales such as the gratings’ period, 
the separation distance or the characteristic wavelength 
of the material m ■ 

Gasimir torques have also been studied between 
anisotropic/birefringent or meta-materials |28H32] where 
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the anisotropy or the misalignement between the two op¬ 
tical axes of two plates in vacuum gives rise to a vacuum 
fluctuation induced torque. The effect is however much 
smaller and we will discuss the orders of magnitude below 
in this paper. 

Casimir torque measurements between nanostructured 
plates could constitute an efficient alternative way 
for testing non-trivial geometrical dependences of the 
Casimir effect. As the lateral force, the torque has the 
advantage that it can be measured without any refer¬ 
ence to the normal Casimir force. However, if the lateral 
Casimir force was measured by a clever adaptation of 
a standard atomic force microscope (AFM) experiment 
[5], torsional motions can not be easily accessed on such 
systems. The experimental strategy must therefore be 
completely different from the one adopted for the lateral 
force. We show later-on that torsion balance techniques 
are perfectly suited to such experimental aims. 

Theoretical expressions for the Casimir energy Ecas 
between corrugated plates have been given, based on a 
general non-specular scattering approach of the Casimir 
effect [T^. These expressions, specifying the angular mis¬ 
match dependence of Ecas, have allowed evaluating in¬ 
duced torques up to second order in the corrugation am¬ 
plitudes m- This result contains the pure lateral force 
as the 6 — 0 special case and the dependence on the mate¬ 
rial properties and corrugation period Ac of both lateral 
forces and torques are determined by identical response 
functions in the perturbative regime m- 

Casimir torques have not yet been computed beyond 
the perturbative regime. But as clearly demonstrated 
theoretically m and experimentally |22[ I23| with nor¬ 
mal and lateral Casimir forces, shape effects turn out to 
be more dramatic when the aspect ratio of the corru¬ 
gations increases. By analogy, new dynamical regimes 
are expected at the level of torques that are induced be¬ 
tween deep corrugations with respect to torques induced 
between shallow ones. This could be particularly true 
when considering metallic plates. It has indeed been re¬ 
alized that the transition from shallow to deep metal¬ 
lic gratings is accompanied by important optical changes 
|34| . Starting with weakly bounded surface plasmon res¬ 
onances defined on shallow gratings, propagating along 
the corrugations, new families of resonances are accessi¬ 
ble on deep gratings which are related to the localization 
of the electromagnetic held inside the grooves [5^. In 
the limit of one-dimensional slit arrays with sub-micron 
dimensions, scattering amplitudes are strongly enhanced 
at resonant wavelengths |36| . A Casimir torque will nec¬ 
essarily be modified by the contributions of such surface 
and guided resonances, suggesting important possibili¬ 
ties for tailoring the dynamical signatures associated to 
Casimir torques. 

We will now present calculations of the Casimir torque 
between two rectangular gratings. Notations of the pa¬ 
rameters are indicated in the schematic view of a single 
grating on Fig. For simplicity the two gratings are 
supposed to have the same geometric parameters. 



FIG. 2. Schematic representation of a single grating with the 
geometrical parameters: a and w are respectively the corru¬ 
gation amplitude and width, d is the grating’s period. 


The Casimir torque r is defined as 

T = —deEcas ( 1 ) 

where Ecas is the Casimir energy and 6 the torsion an¬ 
gle between the two gratings. At zero temperature the 
Casimir energy between the two corrugated plates reads 

m 


Ecas = d^klndet [1 - RiTRaTI . 

( 2 ) 

It is written in terms of the purely imaginary spatial fre¬ 
quency after a Wick rotation. With respect to m we 
have performed a change of variable ^/c —^ here. The 
operator 7~ represents a translation of length L from the 
first to the second grating and 1 is the identity opera¬ 
tor. Ri and R 2 are the reflection operators of the two 
gratings the derivation of which has been detailed else¬ 
where m- The k-integration runs over the first Brillouin 
zone (FBZ). The FBZ is a hexagon in the {kx, ky) recip¬ 
rocal space as depicted on Fig. which is defined as 


<( 


TT , 7r\ 

, 0 27r 0 , , 9 2tt 6 

kx tan —I tan - < < kx tan-tan - 
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(3) 

n 


—kx cot ^ ^ CSC 9 < ky < —kx cot 9 — 



If we choose the lines of the gratings to be along 
the y axis, the first grating will couple via diffraction 
the wavevectors k -|- n’^Gx whereas the rotated grat¬ 
ing will couple the wavevectors k -|- where = 

cos9ex + sin0ej,. The index n belongs to the interval 
[—N, Al] with N the highest diffraction order retained in 
the calculation. 

The common basis between Ri and R 2 , including two 
states of orthogonal polarizations a = {e, h}, is the gener¬ 
alized outer product (k-l-n^ej,)® (k + m^e^)® (e,/i) of 
dimension 2{2N + 1)^ which we will refer to as |n, m, a). 


















This is the basis in which the operators appearing in 
eq. © are represented. For instance, we therefore have 
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{n,m,a\'T\rL ,m ^a') = e 

i^lm = + k: 

27T 


^nn' ^mm' 

2 

nm 


(4) 


= k + — (nea; + me'j.). 




FIG. 3. Schematic structure of the operator Ri in a minimal 
basis |n, m,CT) with = 1 whose size is 2{2N + 1)^ = 18. 
The zeros are represented in white. The lower part shows the 
wavevectors which are coupled by diffraction with red arrows. 
The FBZ is indicated as the grey region. 

For illustration purposes we show in Figs. and 
the schematic structure of the operators Ri and R 2 
when represented in a minimal basis |n,TO,cr) for the 
first diffraction order iV = 1. The first grating rep¬ 
resented by the operator Ri diffracts in the direction 
e^; as shown in the lower part of Fig. Therefore 
(n, m, tT|Ri|n', to', a') oc Smm'- On the other hand, 
the rotated grating represented by the operator R 2 in 
fig-i diffracts only in the direction and we have 
(n, TO, cr|R 2 |n', to', cr') oc Snn'- This leads to the struc¬ 
ture shown in the upper part of Fig. 

The area A of the FBZ is sin 0. The area cov¬ 

ered in k-space by the integration in eq. ([^ is therefore 
(27V + 1)^A. A prohibitively large number of diffraction 
order N becomes necessary at a small torsional angle 
0 in order to sample the values of the integrand which 
contribute to the Casimir energy as A —1 0. This stems 
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FIG. 4. Same as fig. [^for the operator R 2 . 


from the fact that for the combined system of the two ro¬ 
tated gratings, there is a divergent effective period 
in the y-direction. For this reason, small torsional angles 
are numerically unreachable by eq. ([^. We choose to 
restrict our calculations for 6 > 5°. Note that this for¬ 
malism does not gracefully tend towards the formalism 
for unrotated gratings as 0 —>■ 0 since the FBZ does not 
tend towards the FBZq for unrotated gratings which is 
simply ^ < kx < ^,—00 < ky < 00 }. This is due 
to the fact that we consider here gratings with infinite 
dimensions in the x and y directions. 

We have calculated the Casimir energy E{9) as a func¬ 
tion of the torsional angle between two gold gratings. 
The parameters of the gratings are a period d = 400 nm, 
a filling factor / = 0.5 and the height of the corruga¬ 
tions is a = 200 nm leading to an aspect ratio of 1. The 
separation distance is L = 100 nm. The results are 
shown in Figs. © and ©. From these figures, we infer 
a maximum torque per unit area of 3.5 x 10“® N.m“^ 
for a torsional angle 6 « 12.5°. Considering a grating 
surface of ca. 1 mm^, such a level of torque intensity 
appears well within the range of detectability of the tor¬ 
sion balance setup we present below. This comparably 
high value should be contrasted with values of Casimir 
torques induced between anisotropic dielectric plates |28F 
132] which are about 3 orders of magnitude smaller when 
compared at the same distance. We have also per¬ 
formed calculations for larger corrugations height up to 
0 = 2 yia merely leading, somewhat counter-intuitively, 
to a marginal increase in the Casimir torque. 

So far the torque per unit area values reported between 





























































































































































4 



FIG. 5. The Casimir energy Ecas{d) as a function of the 
torsional angle 6 between two gold gratings. 



FIG. 6. The Gasimir torque t{6) as a function of the torsional 
angle 9 between two gold gratings. 

birefringent/anistotropic plates are < 4 x 10“^° N.m“^ 
at distances of the order of 100 nm and thus turn out 
to be too faint for an experimental demonstration of this 
birefringent effect. In a recent paper |32| an experiment 
at much shorter distances (10-50 nm) has been proposed, 
however leading only at the shortest distance of 10 nm 
to a birefringent torque comparable to the one between 
nanostructuerd plates at L = 100 nm. 

In this context, we propose the design of a torsion bal¬ 
ance aimed at measuring weak torques and more specif¬ 
ically Casimir torques between nanostructured plates. 
Our marionette setup is schematically described in Figj^ 
It consists in a small cylindrical plate, the marionette, 
suspended above a grounded plate, the substrate, whose 
relative positioning is controlled via a piezo-stage. With 
nanostructures imprinted on their facing sides, the plates 
will essentially interact through the Casimir energy at 
sufficiently small gaps. Different types of material can be 
considered, from intrinsic Si to metal coated surfaces and 
nanostructures can be milled using for e.g. high resolu¬ 
tion focused ion beam (FIB) or electronic beam lithogra¬ 
phies. 

The marionette in a given plane A is suspended from a 



FIG. 7. The marionette is a suspended thin cylinder, of typi¬ 
cal mass 0.5 g and radius p = 5 mm. The cylinder can be a di¬ 
electric (e.g. intrinsic Si) as such, or coated by a thick layer of 
metal (e.g. Au). One-dimensional corrugations are imprinted 
on the lower side of the marionette which is then suspended 
by ^ = 20 cm from a rigid frame by three threads with a 3-fold 
rotational symmetry (symbolic triangles on the left-hand side 
of the frame). On plane A, the supporting points are sepa¬ 
rated from the symmetry axis of the setup by a and b for the 
hanging points taken on plane B. Typical values can be cho¬ 
sen as a = 1 mm and 6 = 60 ^m. The substrate is an identical 
plate with the same corrugations than that of the marionette. 
Relative spatial and angular positions between the plates are 
controlled by a high-precision piezo-electric stage. The whole 
frame rests on an efficient vibration isolation platform and is 
put in highly quiet environmental conditions. 


fixed frame in plane B by three threads displayed around 
the axis of the balance with a 2tt/3 rotational symme¬ 
try. This tripod geometry is the key feature of our setup 
since it makes it sensitive to torsional degrees of freedom 
while being gravity dominated. Indeed, small rotations 
of our suspended marionette are essentially translated via 
the tripod geometry into gravity pendulum-like displace¬ 
ments of each thread. As a consequence, the potential of 
the suspension is determined by the vertical position of 
the center of mass of the marionette with the important 
consequence that the intrinsic torsional spring constants 
of each thread do not enter into the description of the 
whole suspension setup. This clearly contrasts with clas¬ 
sical torsion balance setups where the potential is actu¬ 
ally governed by the intrinsic torsional spring constant of 
the hanging fiber, putting the effort there on seeking for 
the smallest torsional damping factors and for the largest 
leverages |37l |38] . In our case, the leverage corresponds 
to the lateral extensions of the corrugations themselves, 
lengths that can only be of the order of the mm range. 
This stresses the crucial need for an original and highly 
sensitive design. 

Yet, as for any suspended system, both rotational and 
translational degrees of freedom have to be considered 
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here. Induced forces F and torques r will be inferred 
from lateral displacements ^ and deflection angles 0. Cor¬ 
responding sensitivities will be measured respectively by 
an effective stiffness k = F/^ and an effective torsion co¬ 
efficient K = t/9. In the static case, one has k = Mujf 
and K = loJ^, where M is the mass and I = Mp^/2 the 
moment of inertia of the marionette (a cylinder of radius 
p), LOt and ujr the frequencies of the translation and rota¬ 
tion proper modes. It is important for torque detection 
to favor sensitivity to rotation over translation degrees 
of freedom by decoupling them from each other. This is 
particularly critical in the context of Casimir torques as 
shown in m where the Casimir interaction between two 
corrugated metallic plates can display a complex surface 
potential corresponding to possible mixing between rota¬ 
tions and lateral displacements. This has to be controlled 
in the experiment, and the chosen dimensions of the sus¬ 
pension can already help in this direction. The mass 
and dimensions (M, p) of the marionette, the length i of 
the threads and the positions of the fixation points rel¬ 
ative to the rotational axis of the suspension (a and b 
respectively in plane A and plane B) can indeed be cho¬ 
sen carefully so that ujt ^ oJr- Calculating the geometry 
gives = 2ab/p^, with simply given by g|^, g be¬ 

ing the acceleration of gravity. The decoupling condition 
then corresponds to a,b << p. 

When considering realistic values of the parameters, as 
detailed in the caption of Figj^ we reach a ratio between 
the two frequencies larger than 10 with uJrl2,TT ^ 0.08Hz. 
Having the lowest frequency, the rotation mode is natu¬ 
rally the most easily excited one, and the choice of the 
a parameter gives an additional freedom to tune uJr to a 
specific frequency range. From an experimental point of 
view, it is particularly interesting to be able to set the 
rotational frequency to the lower part of the spectrum, 
below the microseismic fundamental mode at ^ 0.1 Hz, 
hence isolating rotations from most of the seismic back¬ 
ground. 

Provided that by design, the rotational degrees of free¬ 
dom are kept well separated from the translational ones, 
we model our balance as a damped torsional oscillator 
with an effective damping factor 7 . The equation of mo¬ 
tion of the suspended mirror under the influence of both 
thermal fluctuations and the applied Casimir torque is 
thus 


-I-70-I-k 6 > = Tth + TCas- (5) 


Moment of inertia I and torsion coefficient k are de¬ 
termined from the setup dimensions and set the rotation 
proper mode frequency Ur- We can estimate the perfor¬ 
mances of this ideal system starting from noise limita¬ 
tions in the free pendulum. In this case, with no exter¬ 
nal torque applied, the system is only driven by thermal 
fluctuations, and the brownian angular noise spectrum is 
given as 


Se{uj) 


{k — hjj'^Y + ’ 


( 6 ) 


with fee the Boltzmann constant and T the absolute envi¬ 
ronmental temperature. The Casimir torque is measured 
through the angular displacement 0 of the suspended mir¬ 
ror. The estimation of the minimal detectable torque de¬ 
pends on the measurement method and different schemes 
can be considered |39]. We focus here on measuring the 
applied torque by detecting quasi-static changes in the 
angular position of the suspended mirror. The least de¬ 
tectable torque in this situation is given from the low- 
frequency limit of Eq. § as |4(i] 

<5'rmin(w ^ uJr) ~ \/‘2-k^T^ N.m/-\/Hz (7) 

At first sight, the sensitivity of the balance improves as 
the measurement time is increased and the damping fac¬ 
tor reduced. However, achieving better sensitivity by re¬ 
ducing 7 (i.e. increasing the quality factor of the torsion 
balance) can turn detrimental at low frequencies since 
the dynamical response of the now high quality factor 
oscillator will not be stable over such long measurement 
times. In this situation, it becomes necessary to increase 
the bandwidth of the system, through an active control 
of the damping factor. Such an active control (known 
as cold-damping technique) increases the damping factor 
without adding new fluctuations to the system and al¬ 
lows stabilizing the mechanical response of the whole bal¬ 
ance. This technique is efficiently implemented in high- 
precision opto-mechanical |4Hl43] and electro-mechanical 
experiments |44fl46| . 

We can still give simple estimations regarding the po¬ 
tential sensitivity of our setup, assuming, possibly under 
active control, critical damping with 7/2 = hor, in which 
case ^Tmin y/Ak'sTIujr N.m/-\/Hz. With the figures 
given above, this corresponds to 10“^^ N.m/-\/Hz for the 
least detectable torque. This value nicely falls within the 
theoretical evaluation of the Casimir torque given above. 
With actual numbers that will eventually be chosen ex¬ 
perimentally in order to maximize the sensitivity (while 
keeping systematic effects at a minimal level), this pre¬ 
liminary value is already comparable to the sensitivity of 
existing high precision torsion balances EHlIiZl. 

An angular change of the suspended mirror can be 
recorded by an associated displacement of a weak laser 
beam reflected from the center of the mirror onto a low- 
noise position sensitive detector. Such a technique can 
reach a typical sensitivity of a few nrad/ -s/Hz [IS] . This 
detection scheme can be surpassed by optical interfer¬ 
ometers m- Interestingly, interferometers can be made 
sensitive to both rotational and lateral displacements of 
the marionette. Being additionally well adapted for im¬ 
plementing an active damping of the system, the inter¬ 
ferometers can provide a full control over the excitation 
modes of the balance, as they simultaneously allow for 
monitoring the displacements of the marionette and for 
injecting a back-action torque for active control. Also, 
as well-known when operating MEMS and NEMS and 
well-recognized in all Casimir force measurements, bolo- 
metric and electrostatic backgrounds can easily dominate 
to screen the dynamical signal m- They will have to be 
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kept under careful control, though residual spurious ef¬ 
fects might easily be rejected from the data using the 
TT-periodicity of the torque signal m- 

With such a high sensitivity, a first Casimir torque 
demonstration therefore appears within experimental 
reach. With sufficient confidence in the sensitivity, it 
might even become possible to work with relatively large 
plate separations. While this would obviously reduce the 
torque intensity, it would have important experimental 


advantages, gaining over two major technical issues: a 
large working distance reduces tilt effects on the torque 
due to deviation from strict parallelism and minimizes 
the problem of interstitial dust particles, too. Finally, we 
note that here we have limited the discussion to quasi¬ 
static measurements. One has to keep in mind that dy¬ 
namic mode detection can potentially yield better sensi¬ 
tivity, in particular as far as 1// noise limitation is con¬ 
cerned. 
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